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ABSTRACT 


Species of the Chinese endemic genera Yinshania, Hilliella, and Cochleariella were originally placed in or closely 
associated with Cochlearia. A previous preliminary molecular study mainly on European Cochlearia and detailed 
morphological studies by us showed that this complex was not affined to Cochlearia s. str. Depending on the authority 
consulted, the number of taxa recognized in this complex ranged from 11 to 25 species in one to four genera. The 
present phylogenetic study is based on the analysis of the ITS (internal transcribed spacer regions of the nuclear 
ribosomal DNA) and the chloroplast trnL-intron sequences from 18 taxa. Resulting phylogenies were compared, and 
the results demonstrate that there are two different lineages. One lineage combines exclusively the highly polyploid 
taxa from Hilliella and Cochleariella. The second lineage includes the diploid taxa from Yinshania. However, incon- 
gruencies when nrDNA- and cpDNA-derived phylogenies were compared suggest hybridization between these two 
lineages. We followed a concept to combine all taxa of this complex into the genus Yinshania. Our results from 
phylogenetic analysis of nr and cpDNA support the association of Yinshania with Cardamine and Rorippa, rather than 


with Cochlearia, as was suggested by nearly all previous authors. 
Key words: Brassicaceae, Cochleariella, Hilliella, molecular systematics, reticulate evolution, Yinshania. 





Many authors follow Schulz (1936) and Schultze- 
Motel (1986) in dividing Cochlearia into the sec- 
tions Pseudosempervivum Boiss., Glaucocochleria O. 
Ε. Schulz, Cochlearia (= Eucochlearia Prantl), and 
Hilliella O. E. Schulz. As shown by Koch et al. 
(1999a), however, this sectional classification is 
highly artificial. Section Cochlearia is widely dis- 
tributed in Europe and the circumpolar region, 
whereas section Glaucocochlearia, which was 
raised to the generic rank by Pobedimova (1968), 
is restricted to southwestern Europe. The latter sec- 
tion consists of C. glastifolia L. and C. megalos- 
perma (Maire) Vogt, as well as C. aragonensis Coste 
& Soulié, which was only recently included (Koch 
et al., 1996), although considered to be distantly 
related to the other two species (Koch et al., 
1999a). Section Cochlearia consists of a species 
complex that demonstrates highly polymorphic 
chromosome numbers, and diverse ecological ad- 
aptation and geographic distributions. Morphologi- 
cal differences between phylogenetically sister taxa 
are often weak and poorly defined (Koch et al., 
1996). Both sections Cochlearia and Glaucococh- 
learia are closely related to the genus Ionopsidium 


Rchb. (Koch et al., 1999a). Section Pseudosemper- 
vivum, which is centered іп the Middle East and 
clearly unrelated to Cochlearia, is most closely re- 
lated to Masmenia F. K. Mey. and Noccaea Moench 
(Koch et al., 1999a), both of which were segregated 
by Meyer (1973, 1979, 1991) from Thlaspi L. s.l. 
The family of Brassicaceae is divided into ѕеу- 
eral tribes and subtribes. Most of them are highly 
artificial, such as tribe Arabideae (Koch et al., 
1999b) or Lepidieae (Zunk et al., 1996). Following 
classical tribal concepts, Cochlearia sect. Pseudo- 
sempervivum, sect. Cochlearia, and sect. Glauco- 
cochlearia are members of tribe Lepidieae. Species 
originally assigned by Schulz (1923) to section Hil- 
liella (Yinshania, Hilliella, and Cochleariella) were 
excluded from Cochlearia by Pobedimova (1970), 
who did not assign them to any genus. However, 
these species have recently been placed in three 
Chinese endemic genera, Yinshania Y. C. Ma & Y. 
Z. Zhao, Cochleariella Y. H. Zhang & Vogt, and 
Hilliella (O. E. Schulz) Y. H. Zhang, each of which 
was assigned to a different subtribe. The genus Yin- 
shania (Ma & Zhao, 1979) was placed in subtribe 
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Descurainiinae, tribe Sisymbrieae. The genus 
Cochleariopsis (Zhang, 1985), renamed as Coch- 
leariella (Zhang & Cai, 1989), was placed in sub- 
tribe Cochleariinae, tribe Lepidieae, along with Hil- 
liella s. str. (Zhang, 1986). Although a few studies 
on taxonomy, evolution, and origin of these genera 
(Zhang, 1987; Zhang & Xu, 1990; Zhang, 1993) 
have been made, nothing was said about their sys- 
tematic position in relation to the remaining Asian 
and European taxa of Cochlearia. Zhang’s (1987) 
division of Yinshania (excluding Hilliella) into two 
sections and two series and Zhao’s (1992) classifi- 
cation of Yinshania (including Hilliella) into two 
sections and six series, were shown by Al-Shehbaz 
et al. (1998) to be highly artificial. In fact, one of 
the species assigned by Zhang (1987) to Hilliella 
and by Zhao (1992) to Yinshania was placed by Al- 
Shehbaz and Yang (1998) in the synonymy of Car- 
damine fragariifolia O. E. Schulz. 

On the basis of a comprehensive morphological 
survey of Yinshania, Hilliella, and Cochleariella, 
Al-Shehbaz et al. (1998) reduced the latter two to 
synonymy of Yinshania, and concluded that there 
is no need for infrageneric subdivisions that do not 
reflect the phylogenetic relationships of this small 
genus of 13 species. 

They also demonstrated that morphological char- 
acters previously used in the delimitation of species 
(e.g., density of papillae on the fruit valves, fruit 
shape, and seed number per fruit) are highly vari- 
able among and within different populations of the 
same species. Furthermore, differences in the com- 
pression of fruit (terete vs. latiseptate or angusti- 
septate) that were emphasized heavily by earlier 
authors (e.g., Schulz, 1936) in the delineation of 
genera were not found to be taxonomically useful 
in the Yinshania complex. As shown by Koch et al. 
(1999a), the placement of heavy emphasis on fruit 
compression has led to the artificial integration of 
several taxa into Cochlearia sect. Pseudosempervi- 
vum instead of Thlaspi s.l. In fact, terete and var- 
iously flattened fruits occur in numerous genera of 
the Brassicaceae, and in many cases this aspect of 
fruit morphology is taxonomically insignificant. 

Morphological convergence and parallelism are 
widespread in the Brassicaceae (Dvorák, 1971; 
Meyer, 1973; Avetesian, 1983; Endress, 1992), and 
the dependence on such characters to construct 
phylogenies often leads to erroneous conclusions 
(Sytsma, 1990; Meyer, 1991). Recent molecular 
analyses (e.g., Warwick et al., 1992; Price et al., 
1994: Mummenhoff & Koch, 1994; Zunk et al., 
1996; Mummenhoff et al., 1997; Koch et al., 1998a, 
b; Koch et al., 1999a, b) have made significant con- 
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tributions to a better understanding of the classifi- 
cation, generic delimitation, and phylogenetic re- 
lationships in the Brassicaceae. A preliminary 
study (Koch et al., 1999a) utilizing ITS nrDNA and 
cp trnL intron sequence data of four species of the 
Yinshania complex (including Hilliella and Coch- 
leariella) clearly showed that the complex is unre- 
lated to Cochlearia. In this analysis it has been 
shown that ITS апа trnL intron sequence data pro- 
vide sufficient sequence variation to distinguish 
significantly between Yinshania and Cochleariella/ 
Hilliella accessions with both data sets. In order to 
gain a better insight of the phylogenetic relation- 
ships within this Chinese complex, we examined 
sequence variation of the internal spacer regions 
(ITS1 and ITS2) of nrDNA (Baldwin et al., 1995; 
Campbell et al., 1995) and of the cp trnL intron 
(ВӧЫе et al., 1994; Gielly & Taberlet, 1994; van 
Ham et al., 1994; Koch et al., 1999a), and com- 
pared the derived molecular phylogenies with tra- 
ditional concepts based on morphological data. 
This approach provided us with the opportunity to 
characterize species lineages and to analyze incon- 
gruencies between different data sets in order to 
test hypotheses of gene flow over lineages and chlo- 
roplast capture. 


MATERIALS AND METHODS 
PLANT MATERIAL 


Leaf material for DNA extraction was obtained 
from herbarium specimens (Table 1), most of which 
were provided and determined by Zhang Yu-hua 
(Institute of Materia Medica, Zhejiang Academy of 
Medicine, Hangzhou, People’s Republic of China). 
We did not examine vouchers to verify Zhang’s de- 
terminations. The samples represent a broad spec- 
trum of species of Yinshania, Hilliella, and Coch- 
leariella. Cardamine flexuosa With. and Rorippa 
palustris (L.) Besser served as the outgroups. The 
DNA sequences for the outgroups were obtained 
from Franzke et al. (1998). 


DNA EXTRACTION, PCR-AMPLIFICATION, AND 
SEQUENCING 


The total DNA for the outgroups was isolated 
from leaf tissues following the CTAB (cethyltriam- 
moniumbromide) method of Doyle and Doyle 
(1987), as modified by Mummenhoff and Koch 
(1994). DNA extraction from herbarium material 
was performed in a mini preparation in Eppendorf 


reaction tubes from 50-ywg dried tissue. Tissue was 
ground with sand and prewarmed 2X CTAB-buffer. 
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Organic extraction and DNA isolation from herbar- 
lum specimens followed Koch et al. (1996). 

Double-stranded DNA of the complete ITS re- 
gion, including the 5.8S rDNA gene, was amplified 
by 30 cycles of symmetric PCR using ITS primers 
Initially designed by White et al. (1990) and mod- 
ified by Mummenhoff et al. (1997). The 18F primer 
(5'-ССААССАСААСТССТААСААСС-3') is lo- 
cated at the 3'-end of the 18S rDNA gene, and 
primer 25R (5'-TCCTCCGCTTATTGATATGC-3') 
is located at the 5'-end of the 25S rDNA. It has 
been reported that PCR selection of rDNA para- 
logues has occurred (Buckler et al., 1997). How- 
ever, PCR selection might have only been impor- 
tant in high G+C content sequences (Buckler et 
al., 1997). Sequences from Yinshania and Hilliella 
(Koch et al., 1999a) are comparable in G+C con- 
tent to sequences from Gossypium, in which PCR 
selection was probably weak (Wendel et al., 1995a). 
The resulting amplification product included ITS1, 
5.8S rDNA, and ITS2. Only those PCR products 
were cloned into the pGEM-T-Easy cloning vector 
(PROMEGA) that showed a single band on ethi- 
dium bromide stained agarose gels. Two cloned ITS 
regions from two independent PCR reactions were 
sequenced (forward and reverse) with both ampli- 
fication primers and two universal primers located 
in the flanking sites of the pGEM-T-Easy vector (t7- 
forward: 5'gtaacgatttaggtgacactatcg-3, m13-re- 
verse: 5'-agcggataacaatttcacacagga-3). This means 
that every single clone was sequenced four times 
to avoid sequence errors. The trnL (ОАА) intron 
was amplified and sequenced by using the univer- 
sal primer B49318 (5'-ССАААТСССТАСАСССТ- 
ACG-3') located at the 3'-end of the trnL(UAA)5'- 
exon and А49855 (5'-ССССАТАСАСССАСТТС- 
AAC-3") located at the 5’-end of the trnL(UAA)3’- 
exon (Taberlet et al., 1991). The PCR profile used 
to amplify the trnL intron followed the following 
profile: hot start with 5 min. at 94°C, and 35 cycles 
of amplification (1 min. 94°C, 45 min. 50°C, 45 
min. 72°C), final elongation step for 10 min. 72°C, 
and storage at 4°C. DNAs were cycle-sequenced 
using the Taq DyeDeoxy Terminator Cycle Se- 
quencing Kit (ABI Applied Biosystems, Inc.). Prod- 
ucts of the cycle sequencing reactions were run on 
an ABI 377XL automated sequencer (ABI Applied 
Biosystems, Inc.). Material from accession numbers 
22-26 (Table 1) was only used for sequencing the 
trnL intron sequence, because amplification of the 
ITS regions failed totally. 


DATA ANALYSIS 


ITS data. Boundaries of the ITS regions and 
coding sequences were determined by comparison 
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to those of Sinapis alba L. (Rathgeber & Capesius, 
1989) and other Brassicaceae (Mummenhoff et al., 
1997; Koch εἰ αἰ., 1999а). DNA sequences were 
aligned by hand. Parsimony analyses меге рег- 
formed with unordered Fitch parsimony and 
weighted parsimony with a transition: transversion 
weighting of 1.0:1.08 using PAUP version 3.1 
(Swofford, 1993). The BRANCH-AND-BOUND al- 
gorithm was used to find maximally parsimonious 
trees. Bootstrap analysis (Felsenstein, 1985) was 
performed using 1000 replicates and the HEURIS- 
TIC search algorithm with the MULPARS option. 
We combined the GAPMODE=MISSING option 
with the coding of the gaps as additional presence/ 
absence characters (Downie & Katz-Downie, 1996). 
This option decreases the number of equally par- 
simonious trees because of the redundancy result- 
ing from having two sets of scored characters for 
the same indel events (Wojciechowski et al., 1993). 
Evolutionary data are most often presented as a 
phylogenetic tree, the underlying assumption being 
that evolution is a branching process. However, em- 
pirical data is rarely ideal and often supports sev- 
eral trees instead of one unique tree. Hence, it 
makes sense to consider tree reconstruction meth- 
ods that produce a tree if the given data heavily 
favor one tree over all others. Otherwise, methods 
that produce a more general graph that indicates 
different possible phylogenies are useful (Huson, 
1998). One such method is the Split Decomposition 
introduced by Bandelt and Dress (1992) and its 
variations. In order to visualize conflicting phylo- 
genetic signals, we analyzed all ingroup ITS se- 
quences (see Fig. 6), using the software program 
SplitsTree version 1.0.3. (Huson & Wetzel, 1995 
shareware ftp://ftp.uni-bielefeld.de/pub/math/splits/). 
trnL intron data. DNA sequences were aligned 
by hand. Parsimony analysis was performed with 
PAUP (version 3.1; for options, see ITS data). Gaps 
were treated as additional unweighted binary char- 
acters. These gaps were coded using strict criteria: 
gaps must occur at the same position and have the 
same aligned length to be treated as homologous, 
and no splitting of one gap in two or more char- 
acters was performed (Koch et al., 1999a). Boot- 
strap analysis was performed as described above. 


TRIBAL RELATIONSHIPS 


To estimate the tribal affinity of Yinshania, we 
derived an ITS phylogeny with sequences from 
Capsella rubella Reut. (Koch et al., 1999b), Ara- 
bidopsis thaliana (L.) Heynh. (GenBank U43224), 
Yinshania acutangula (O. E. Schulz) Y. H. Zhang, 
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Barbarea vulgaris R. Br. (EMBO X98632), Carda- 
mine flexuosa With. (Franzke et al., 1998), Thlaspi 
arvense L. (Koch et al., 1999a), Cochlearia aestu- 
aria (Lloyd) Heywood (Koch et al., 1999a), Brassica 
oleracea L. (GenBank AF039994/AF040038), and 
Sinapis alba (EMBO X66325). DNA sequences 
were aligned by hand, and the alignment is shown 
in Figure 4. Alignment positions 109-160 were re- 
moved from subsequent analysis. Parsimony anal- 
yses were performed with unordered Fitch parsi- 
mony using PAUP version 3.1 (Swofford, 1993). 
The BRANCH-AND-BOUND algorithm was used 
to find maximally parsimonious trees with the GAP- 
MODE=NEWSTATE option. Bootstrap analysis 
(Felsenstein, 1985) was performed using 1000 rep- 
licates and the HEURISTIC search algorithm em- 
ployed with the MULPARS option. A decay anal- 
ysis (Bremer, 1988) was performed in addition to 
the bootstrap approach, in order to assess the con- 
fidence that could be placed in the monophyly of 
clades. Decay indices (DI) were estimated accord- 


ing to Baum et al. (1994). 


RESULTS 
ITS DATA 


The total length of the alignment with accession 
numbers 1—21 and outgroups Cardamine flexuosa 
and Rorippa palustris is 464 bp, with 283 and 181 
nucleotides in ITS] апа ITS2 spacer regions, re- 
spectively. The alignment required 36 (7.8%) gap 
positions, including the outgroups, and is shown in 
Figure 1. Sequence data were submitted to Gen- 
Bank with accession numbers ΑΕΙΟΟΤΟΣ-- 
AF 100852 (Table 1). One third of these gaps is 
located between positions 117 and 133 in the ITS1. 
This region was completely excluded from subse- 
quent analysis. Additionally, we excluded nucleo- 
tide position 465-467 and 473-484 from subse- 
quent data analysis because of an ambiguous 
alignment (Fig. 1). Gaps from nucleotide positions 
11-14 and 293-294 were treated as one single gap, 
respectively. The number of introduced gaps 15 
comparable to a phylogenetic analysis of Thlaspi 
s.l. with 4.8% of the sites (Mummenhoff et al., 
1997). In Krigia and outgroups, Kim and Jansen 
(1994) had to introduce gaps in 3.9% of the sites. 
Total lengths of [Τ5] and ITS2 are nearly identical 
among the taxa surveyed, and vary between 447 bp 
(Hilliella alatipes (Hand.-Mazz.) Y. H. Zhang & H. 
W. Li var. micrantha Y. H. Zhang [= Yinshania 
rivulorum (Dunn) Al-Shehbaz et al.], accession no. 
10) and 457 bp (Cardamine flexuosa and Rorippa 
palustris). 

Phylogenetic analysis using Fitch parsimony, in- 
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cluding the additional 0/1 matrix for the gap posi- 
tion, resulted in 24 most parsimonious trees (MPTs) 
with a length of 538 and a consistency index (CI) 
of 71.4% (66.9% if autapomorphies are excluded). 
Of the 264 variable nucleotide positions, 166 in- 
formative positions were in the ITS] region (in- 
cluding 48 autapomorphies) and 98 in the ITS2 
region (including 25 autapomorphies). Four out of 
36 gap positions within the total sequence align- 
ment are unique to a particular sequence (Rorippa 
palustris 2 gaps, acc. no. 3, and acc. no. 14). A 
calculation of the transition/transversion ratio for 
the МРТ» revealed a ratio of 1.00: 1.08. Therefore, 
we used a weighted parsimony approach with a 
character state weighting of 1.00: 1.08 (transition: 
transversion), which resulted in one MPT that is 
also represented among the 24 MPTs from the Fitch 
parsimony with a consistency index of 67.7% (CI 
57.9% if autapomorphies were excluded). We pre- 
sent the MPT from the weighted parsimony ap- 
proach to demonstrate relative branch length (Fig. 
2). Bootstrap values are provided from 1000 repli- 
cates using the weighted parsimony approach. For 
most taxa we identified only one ITS sequence type 
within a single specimen. For Hilliella lichuanensis 
Y. H. Zhang (acc. no. 13), Cochleariella zhejian- 
gensis (Y. H. Zhang) Y. H. Zhang & R. Vogt (acc. 
no. 21), Yinshania acutangula (acc. no. 5), Y. henryi 
(Oliv.) Y. H. Zhang (acc. no. 7), and Y. furcatopilosa 
(K. C. Kuan) Y. H. Zhang (acc. no. 8), we found 
two very similar ITS sequences among the two 
clones sequenced. In the case of accession num- 
bers 13, 5, 7, and 8, ITS1 and ITS2 regions differed 
by only a single site mutation. In the case of C. 
zhejiangensis (acc. no. 21), the two ITS types from 
the same individual differed by 15 mutations, 
which might indicate that two different ITS loci 
were cloned and sequenced. Both sequences clus- 
tered within the same clade. In Hilliella fumarioi- 
des (Dunn) Y. H. Zhang & H. W. Li (acc. no. 17), 
we detected two ITS sequences from a single in- 
dividual that clustered in different clades (Figs. 2, 
6). Both sequences differed by 76 mutations. The 
ITS data clearly support a separation of two clades 
consisting of Yinshania sensu Zhang (1987) and 
Hilliella/Cochleariella. Different accessions from 
one taxon (sensu Al-Shehbaz et al., 1998, refer to 
fig. 2) grouped at different positions in the case of 
C. zhejiangensis (acc. no. 20, not related to acc. nos. 
18, 19, 21). 

Using the taxonomic concept of Al-Shehbaz et 
al. (1998) and merging Hilliella changhuaensis Y. 
H. Zhang, H. guangdongensis Y. H. Zhang, and H. 
lichuanensis (acc. nos. 11, 12, 13, respectively) into 
H. lichuanensis, only accession numbers 11 and 13 
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grouped together by ITS data. Accession number 
12 (H. guangdongensis) is unrelated to these ac- 
cessions based on ITS data. Similarly, H. warburgii 
(O. E. Schulz) Y. H. Zhang & H. W. Li (асс. nos. 
18, 19), Cochleariella zhejiangensis (acc. nos. 20, 
21), and H. fumarioides (acc. no. 17) merge into Y. 
fumarioides (Dunn) Y. Z. Zhao, as proposed by Al- 
Shehbaz et al. (1998). ITS sequences from Y. fu- 
marioides sensu Al-Shehbaz et al. (1998) are found 
in three different positions among the Hilliella/ 


Cochleariella clade (Fig. 2). 


trnL DATA 


The alignment of 514 bp is interspersed with 
seven gaps as shown in Figure 3. Sequence data 
were submitted to GenBank with accession num- 
bers AF100853 through AF100881 (Table 1). The 
lengths of trnL intron sequences range from 311 bp 
in Hilliella rivulorum to 514 bp in Cardamine flex- 
uosa. Of the 38 variable nucleotide positions in the 
alignment, 12 sites are autapomorphic. In addition, 
one of the seven gaps is autapomorphic. Phyloge- 
netic analysis resulted in two MPTs with a length 
of 53 steps and a consistency index of 92.5% 
(90.0% if autapomorphies were excluded). One 
most parsimonious tree is shown to demonstrate rel- 
ative branch length. The strict consensus tree could 
be generated easily by drawing branches that are 
indicated in Figure 2 by heavy bars with zero 
length. Bootstrap values are given from 1000 rep- 
licates. 

Discrimination of a Yinshania clade from a Coch- 
leariella/Hilliella clade is not as obvious with trnL 
intron data as compared to ITS sequence data (Fig. 
2). In the trnL tree the branch setting of the Yin- 
shania clade as a sister group to the Hilliella/Coch- 
leariella clade is not highly supported, and boot- 
strap value for this branching point is less than 
50% (Fig. 2). Removing additional gap characters 
from the matrix parsimony analysis resulted in one 
MPT with the Yinshania clade as a sister group to 
H. hunanensis (acc. no. 16), H. guandongensis (acc. 
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Figure 2. 
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no. 12), C. zhejiangensis (acc. no. 20), and H. ala- 
tipes var. micrantha (acc. no. 10). Remaining Hil- 
liella/Cochleariella taxa appeared in this analysis 
for trnL data excluding gap information as a sister 
group to these two clades. Nonetheless, integration 
of Yinshania gianningensis Y. H. Zhang into the 
Hilliella/Cochleariella clade is significant for trnL 
data, in contrast to its segregation by ITS. Esti- 
mation of decay indices (DI) using the trnL intron 
matrix without gap information revealed a high val- 
ue, DI = 3+, for the branch setting of the Yin- 
shania clade. Different accessions of C. zhejiangen- 
sis (acc. nos. 20, 21) did not group closely together; 
this has also been documented for the ITS data 
(Figs. 2, 6). Hilliella changhuaensis, H. guangdon- 
gensis, and H. lichuanensis do not group together 
in the cpDNA-based tree as proposed by the mor- 
phology-based concept combining them in H. li- 
chuanensis (Al-Shehbaz et al., 1998); the same dis- 
cordance holds for H. warburgii, C. zhejiangensis, 
and H. fumarioides. Based on the trnL sequence 
data, they are not combined in one single clade that 
could be named as H. fumarioides as proposed in 
the revison of Al-Shehbaz et al. (1998). 


ITS VERSUS TRNL INTRON DATA 


Both phylogenetic trees (ITS vs. trnL data) from 
Figure 2 show some congruencies (all following ar- 
guments are also true when comparing the strict 
consensus trees from Fitch parsimony, which could 
be easily deduced by drawing branches indicated 
by heavy bars with zero length): 

(1) Hilliella guangdongensis [sensu Zhang], 
which has been merged in H. lichuanensis sensu 
Al-Shehbaz et al. (1998), is separated from remain- 
ing Н. lichuanensis sensu Al-Shehbaz et al. (1998), 
and it is more closely related to H. hunanensis, a 
taxon that was recognized by Al-Shehbaz et al. 
(1998). 

(2) Hilliella warburgii and Cochleariella zhejian- 
gensis are not integrated into H. fumarioides sensu 


Al-Shehbaz et al. (1998), but (3) most accessions 


Comparison of the ITS-derived phylogeny with those from the trnL sequence data from Yinshania, Coch- 


leariella, and Hilliella, as well as outgroups Cardamine and Rorippa from subtribe Arabideae (Brassicaceae). Accession 
enumeration in brackets follows Table 1. For the ITS and the trnL tree the nomenclature sensu Zhang (1985, 1986, 
1987, 1993) has been used. In between both trees the taxonomic treatment of these taxa according to Al-Shehbaz et 
al. (1998) is shown. The most parsimonious tree (weighted parsimony approach) from ITS nuclear sequence data and 
one out of two МРТ» from the trnL plastid sequence data is shown to demonstrate relative branch length. Weakly 
supported branches that collapsed in the strict consensus trees using Fitch parsimony are highlighted by a broad line. 
(The strict consensus trees could therefore be easily generated by drawing broad lines as zero branches.) Bootstrap 
values are provided from 1000 replicates and shown above branches if greater than 50%. Distance scale indicating 
branch lengths is given below. ITS designation (a) and (b) of accession nos. 5, 7, 8, 13, 17, and 21 indicate different 


ITS types from a single individual. 
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of Н. warburgii and С. zhejiangensis formed one 
single cluster with no separation of H. warburgii 
versus C. zhejuangensis. This is in agreement with 
previous concepts combining both taxa, C. zhejian- 
gensis and H. warburgii, in one single taxon. These 
congruent ITS and trnL intron findings demonstrate 
that both Y. lichuanensis and Y. fumarioides species 
complexes, treated as well-defined taxa according 
to the morphological revision of Al-Shehbaz et al. 
(1998), do not form monophyletic groups by molec- 
ular evidences. 

However, some incongruencies among the two 
molecular data sets could be detected: (1) Yinshan- 
ia gianningensis grouped either inside the Hilliella/ 
Cochleariella clade (trnL intron data) or into Yin- 
shania s. str. (ITS data); (2) H. rupicola (D. C. 
Zhang & J. Z. Shao) Y. H. Zhang also clustered 
into two different subgroups within Hilliella/Coch- 
leariella. 


TRIBAL RELATIONSHIPS 


The alignment of ITS1-5.8SrDNA-ITS2 is 644bp 
in length (Fig. 4). Within the ITS1 region (bp l- 
301) there are 143 variable nucleotide positions 
(including 62 autapomorphies). The ITS2 region 
(position 456—644) contains 85 variable nucleotide 
positions (including 57 autapomorphies). The 5.85 
rDNA gene, located between both spacer regions 
(bp 302—455), contains 11 variable nucleotide po- 
sitions (including 6 autapomorphies). Because of an 
ambiguous alignment, nucleotide positions 107— 
153 and 460-484 were removed from the original 
data matrix (Fig. 4), resulting in a final data matrix 
of 412 bp. 

Fitch parsimony analysis resulted in one MPT 
with a length of 288 steps and a consistency index 
of 79.5% (66.5% if autapomorphies are excluded). 
The phylogenetic tree (Fig. 5) showed closer rela- 
tionships of Yinshania to taxa from tribe Arabideae 
sensu Janchen (1942) (Arabidopsis (DC.) Heynh., 
Barbarea R. Br., and Cardamine L. were used to 
represent tribe Arabideae). However, Cochlearia 
(including species loosely affined to the Yinshania 
complex) and Capsella rubella were placed by 
Janchen in the tribe Lepidieae, where they had 
been put by Hayek (1911) and Schulz (1936). Our 
ITS data do not support this latter placement. A 
molecular analysis of Arabidopsis, Arabis L., and 
their relatives shows a close relationship of Cap- 
sella rubella to Arabidopsis thaliana (Koch et al., 
1999b). Also demonstrated is the polyphyly of Ar- 
abis and Arabidopsis, indicating that tribal struc- 
tures within Brassicaceae are highly artificial. The 
ITS phylogeny does suggest that Chinese Yinshania 
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and related taxa are closer to Cardamine and Bar- 
barea from tribe Arabideae than to genera like 
Thlaspi or Cochlearia s. str. of the tribe Lepidieae. 


CYTOLOGY 


Little is known about the cytology of Yinshania. 
Zhang (1995, and pers. comm.) counted 2n = 12 
for Y. qianningensis Y. Н. Zhang [= Y. acutangula 
sensu Al-Shehbaz εἰ al., 1998], Y. henryi (Oliv.) Y. 
Н. Zhang, and Y. furcatopilosa (K. C. Kuan) Y. Н. 
Zhang, 2n = 42 for Hilliella yixianensis Y. H. 
Zhang, H. paradoxa (Hance) Y. H. Zhang & H. W. 
Li, and H. changhuaensis [= Y. lichuanensis sensu 
Al-Shehbaz et al., 1998], and 2n = 44 for H. 
shuangpaiensis Z. Y. Li [= Υ rupicola sensu Al- 
Shehbaz et al., 1998]. These data correspond to the 
ITS-derived phylogeny, in which the diploid Y. 
qianningensis, Y. henryi, and Y. furcatopilosa, to- 
gether with Y. acutangula, are separated from the 
polyploid Hilliella. Within polyploid Hilliella, the 
Hilliella taxa with 2n = 42 are combined. Hilliella 
shuangpaiensis (represented in this study by acc. 
no. 4, Fig. 2) with 2n = 44 did not group closely 
to the known 2n = 42 taxa (represented in this 
study by асс. nos. 2, 11, and 15, Fig. 2). However, 
any conclusions based on cytology must be only 
preliminary. It remains possible that the Hilliella/ 
Cochleariella group could also be represented by 
polyploid taxa with 2n = 42 and 44. Base chro- 
mosome number for this hexaploid group is x = 7, 
instead of x = 6 as in the Yinshania group. A few 
taxa within the Hilliella/Cochleariella clade may be 
aneuploids (2n = 44) that derived from 2n = 42. 


DISTRIBUTION 


Some interesting features emerge when topolo- 
gies of the phylogenetic trees are compared to the 
geographic distribution of Yinshania s.l. Geograph- 
ic distribution of ITS sequence types from taxa of 
the four main clades within Hilliella/Cochleariella 
do not follow their phylogenetic relationships (Fig. 
7), and they are randomly mixed in Southeast China 
in the provinces of Guangdong, Jiangxi, Zhejang, 
Anhui, Hubei, and eastern Sichuan. However, they 
are separated geographically and phylogenetically 
from the Yinshania clade from Hunan, Xizang, and 
Sichuan. Taxa from the Yinshania clade extend the 
distribution to the southwest. A mixed distribution 
of DNA types also holds for the trnL data. No geo- 
graphic structuring of plastome types could be ob- 
served among taxa from the Hilliella/Cochleariella 
clade (Fig. 8). Based оп trnL intron data Y. qian- 
ningensis (acc. no. 1) from Sichuan integrates into 


the Hilliella/Cochleariella clade. Geographically, 
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Figure 5. Phylogenetic relationships of cruciferous taxa from different tribes. Genetic distances were calculated 
using the program TREECON (van de Peer & de Wachter. 1994) under the Kimura model (Kimura, 1980). Gaps were 


not taken into account, and the neighbor joining algorithm was used to calculate genetic distances. The robustness of 
the tree was tested by 1000 bootstrap replicates. 
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Figure 6. Split Decomposition of the ITS sequence data set analyzing the ingroup taxa. The DRAW-EQUAL-EDGES 
option was used to draw the network, and thus distances are not drawn to scale. 
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Distribution of ITS types among Yinshania and Hilliella/Cochleariella accessions under study in China. 


Taxon enumeration follows Table 1. The Hilliella-Cochleariella clade is separated into four ITS types indicated by 
circles and boxes (gray or white) well suppported by high bootstrap values in Figure 2. Accessions marked by open 
circles and boxes are two well-separated subgroups, which are combined to a single group in Figure 8 (plastid trnL 
data). Phylogenetic relationships are shown schematically in the upper left box and follow Figure 2. For the Yinshania 
clade phylogenetic relationships are shown to demonstrate the position of Yinshania qianningensis |Y. acutangula sensu 


Al-Shehbaz] (асс. no. 1) (refer to Fig. 2). 


this accession lies close to a similar plastome type 
(acc. no. 15) from the Hilliella/Cochleariella clade 
(Fig. 8). Upon comparison of the ITS and trnL in- 
tron-derived phylogenies, taxa that showed differ- 
ent positions in the phylogenetic analysis (e.g., H. 
rupicola, acc. no. 14; H. fumarioides, acc. no. 17: 
one additional ITS sequence copy) were located at 
the center of the Hilliella/Cochleariella distribu- 
tional area. We interpret these results as a first bio- 
geographical documentation of reticulation within 
polyploids. Hilliella fumarioides (acc. no. 17) pos- 
sessed a trnL intron type similar to samples from 
surrounding areas. The two ITS types found in that 
particular individual were also present in adjacent 
regions and from different species. From Hilliella 
rupicola (acc. no. 14) we could isolate an ITS DNA 
type, which is also present in adjacent regions in 
accession numbers 10, 16, 17, and 20. trnL intron 
DNA type from accession number 14 also corre- 


sponds to trnL intron DNA types from 2, 3, 11, 13, 
17, and 20. 


MORPHOLOGICAL VARIATION AND TAXONOMIC 
CONSIDERATIONS 


As in numerous other cases in the Brassicaceae, 
morphological differentiation among and within so- 
called genera and taxa of the Yinshania s.l. com- 
plex does not provide ample characters to draw 
sharp and uncontroversial generic boundaries. 
Flower morphology in the entire complex is of no 
predictive diagnostic value, even for the separation 
of species (Al-Shehbaz et al., 1998). Only leaves 
and fruits offer characters useful for the separation 
of species, and all taxa appear to be not well de- 
fined. As shown by Al-Shehbaz et al. (1998), seed 
sculpture and number per locule, cotyledonary po- 
sition, development of the fruit septum, type (if any) 
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Figure 8. Distribution of trnL intron types among Yinshania and Hilliella/Cochleariella accessions under study in 
China. Taxon enumeration follows Table 1. The Hilliella-Cochleariella clade is separated into three plastid trnL types 
indicated by circles and boxes (gray or white) well suppported by high bootstrap values in Figure 2. Phylogenetic 
relationships are shown schematically in the upper left box and follow Figure 2. For the Yinshania clade phylogenetic 
relationships are shown to demonstrate the position of Yinshania gianningensis |Y. acutangula sensu Al-Shehbaz] (acc. 


no. 1) (refer to Fig. 2). 


of fruit compression, fruit shape, development of 
papillae on the fruit valve, and trichome type, all 
are unreliable in dividing the complex into the gen- 
era Yinshania, Hilliella, and Cochleariella. In fact, 
no single character or set of characters can be re- 
lied upon to subdivide the complex into generic or 
infrageneric taxa. The recognition of a single genus 
is, therefore, taxonomically expedient, as there is 
no single character setting one clade apart from 
another. Even pustules on the testa surface of the 
valves, which were assumed to be a good character 
to separate Yinshania (pustules present) from Hil- 
liella (Zhang & Xu, 1990), appear not to be con- 
gruent, failing to provide a good argument to split 
Yinshania into several genera (Al-Shehbaz et al., 
1998). 


MORPHOLOGICAL CLUSTER ANALYSIS 


A cluster analysis of morphological and ecolog- 
ical characters has been performed by Zhang and 


Xu (1990). This analysis considered distribution, 
altitude, and habit, but also morphological char- 
acters describing hairs, inflorescences, flower de- 
tails, silicles, seeds, and cotyledons. They conclud- 
ed that minute pustules on the testa surface of the 
valves were a good character to separate Yinshania 
from Hilliella. Comparison of our molecular phy- 
logenies with this cluster analysis based on 32 mor- 
phological and 3 ecological characters demonstrat- 
ed the lack of congruency between morphological 
and molecular evolution. This cluster analysis sep- 
arated Yinshania sensu Ma and Zhao (1979) from 
Hilliella sensu Zhang and Li (Zhang, 1986) (a sim- 
ilar resolution to the ITS data). Taxa such as H. 
paradoxa, H. lichuanensis, and H. changhuanensis 
are closely related to each other (as confirmed by 
molecular data by acc. nos. 11, 13, and 15 herein). 
However, this morphological cluster analysis also 
grouped H. guangdongensis (acc. no. 12 herein) 
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into this group (Y. lichuanensis clade, Fig. 2). For 
both ITS апа trnL intron sequence, Н. guangdon- 
gensis is very divergent from others grouping with 
the Y. lichuanensis clade (Fig. 2). Hilliella alatipes 
var. micrantha and H. sinuata are separated with 
both molecular markers (Fig. 2), but are combined 
as a well-supported group in morphological cluster 
analysis (Zhang & Xu, 1990). In summary, within 
Hilliella there is little morphological agreement 
with our molecular data: (1) H. paradoxa, H. chan- 
ghuaensis, and H. lichuanensis are closely related 
to each other (as Y. lichuanensis clade, Fig. 2); (2) 
Н. hunanensis and Н. rupicola grouped together on 
the ITS tree (but not based on trnL data, Fig. 2). 
No further correlation could be observed. None of 
the three data sets (ITS, trnL intron, morphology) 
is powerful enough to elucidate phylogenetic signal 
for the whole species complex. However, significant 
correlations could be observed when cpDNA- and 
nrDNA-derived phylogenies were compared, divid- 
ing the Hilliella/Cochleariella clade into several 
subgroups (Fig. 2) clearly separating Yinshania rel- 
atives from Hilliella/Cochleariella with the sole ex- 
ception of Y. gianningensis. 


HYBRIDIZATION, INTROGRESSION, CHLOROPLAST 
CAPTURE, AND CONCERTED EVOLUTION 


The phylogeny based on plastid trnL intron se- 
quence data reflects the maternal lineages because 
plastids are inherited maternally in most angio- 
sperms, including the Brassicaceae (Harris & In- 
gram, 1991; Reboud & Zeyl, 1994). Introgression 
of a chloroplast type characteristic for the Hilliella/ 
Cochleariella clade into Yinshania qianningensis 
demonstrates possible gene flow between both 
groups. Because of the highly polyploid genomes of 
Н. yixianensis, Н. paradoxa, Н. changhuaensis, and 
Н. shuangpaiensis (represented by acc. nos. 2, 15, 
11, and 4, respectively, in Fig. 2) with multiple 
rDNA loci, and the assumed hybridization within 
the Hilliella/Cochleariella clade and even with the 
Yinshania clade, there is a high probability of con- 
certed ITS sequence evolution. In principle, there 
are three different ways that two different ITS cop- 
ies evolve within a single individual: (1) unidirec- 
tional concerted evolution leads to the loss of one 
copy and fixation of the second (detected in H. rup- 
icola acc. no. 14, herein; and in Gossypium, Wendel 
et al., 1995a); (2) both ITS copies are still present, 
which might be mostly the case in young hybrido- 
genous taxa (detected in H. fumarioides acc. no. 17, 
and С. zhejiangensis acc. no. 21, Fig. 2; in Krigia, 
Kim & Jansen, 1994; in Arabidopsis, O’Kane et al., 
1996); and (3) concerted evolution leads to a new 


Annals of the 
Missouri Botanical Garden 


ITS type that represents a mixture of the two orig- 
inal ITS sequences (in Gossypium, Wendel et al., 
1995b; in Microseris, van Houten et al., 1993; in 
Microthlaspi, Mummenhoff et al., 1997). The third 
type of concerted evolution might have happened 
in Н. sinuata (acc. no. 3). This accession showed a 
plastome type more similar to H. shuangpaiensis 
(acc. no. 4). However, ITS sequence types from pu- 
tative parental ITS sequence types from H. sinuata 
(асс. no. 9) and H. shuangpaiensis (acc. no. 4) ex- 
hibit some additive features found in H. sinuata 
(acc. no. 3). Comparing these three sequences, 
there are 39 variable nucleotide positions (21 in 
ITS1 region, and 18 in ITS2 region). Within the 
ITS1 region, 17 out of 21 variable nucleotide po- 
sitions (81%) are identical among the two H. sin- 
uata accessions (nos. 3, 9); H. shuangpaiensis 
shared only three mutations with H. sinuata (acc. 
no. 9) and one mutation with H. sinuata (acc. no. 
3). Within the ITS2 region both H. sinuata acces- 
sions have only one nucleotide position out of 18 
variable nucleotide sites invariant, but H. shuang- 
paiensis (acc. no. 4) shared 15 positions (83%) with 
H. sinuata (acc. no. 3) and two positions with H. 
sinuata (acc. no. 9). These findings indicate that 
ITS type of H. sinuata (acc. no. 3) consists of a 
mixture of ITS types from relatives of H. shuang- 
paiensis (mostly ITS2 region) and H. sinuata (ITS1 
region). The overall sequence divergence between 
H. shuangpaiensis (acc. no. 4) and H. sinuata (acc. 
по. 9) is 7.8%. 

Concerted evolution of ITS DNA loci has been 
shown several times to occur in the Brassicaceae 
(O’Kane et al., 1996; Mummenhoff et al., 1997; 
Koch et al., 1998b; Franzke et al., 1998) and other 
families (Wendel et al., 1995a, b; Buckler et αἰ.. 
1997). Sequence divergence values of ITS types 
from putative parents, giving rise to hybrids in 
which concerted evolution has been observed, 
ranged from 3.1% (Microthlaspi natolicum vs. M. 
perfoliatum, Mummenhoff et al., 1997), 5.0% (Car- 
damine атага vs. С. rivularis auct., Franzke et al., 
1998), to 6% (Arabidopsis thaliana vs. A. arenosa, 
O’Kane et al., 1996) comparable to a value of 7.8% 
found among H. shuangpaiensis versus H. sinuata. 

We conducted a split decomposition analysis to 
visualize conflicting phylogenetic signal indicating 
concerted evolution in groups among the Hilliella/ 
Cochleariella clade (Fig. 6). This analysis clearly 
indicates hybridization with subsequent concerted 
evolution of ITS regions in H. sinuata (acc. no. 3). 
Concerted evolution of ITS sequences greatly influ- 
ences any interpretation of the ITS phylogeny. 
Since diploid members of the Brassicaceae such as 
Arabidopsis thaliana typically show 2 NOR loci, 
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one could assume that in hexaploid Hilliella/Coch- 
leariella taxa at least 6 major NOR loci are present. 
Therefore, sequencing of two individual ITS clones 
is not a sufficient survey to find all putative ITS 
types from a single individual. This undersampling 
leads to an underestimation of ITS type variation 
as well as the degree of hybridization and concerted 
evolution. 

We found evidence suggesting concerted evolu- 
tion and described three examples (H. fumarioides 
ace. no. 17, H. rupicola acc. no. 14, H. sinuata acc. 
no. 3) of possible hybridization and subsequent 
concerted evolution of ITS sequence. The overall 
amount of sequence divergence between taxa in 
this study demonstrates a relatively high age for the 
different lineages. In fact, sequence divergence val- 
ues are much higher when compared to those ob- 
tained in infrageneric studies of closely related spe- 
cies of other Brassicaceae (e.g., Cochlearia s. str. 
< 1.75%, Koch et al. (1999a); Noccaea < 4.5%, 
Mummenhoff et al. (1997); Cardamine < 4.5%, 
Franzke et al. (1998)). Within the Yinshania clade, 
sequence distance values range up to 6%, but with- 
in the Hilliella/Cochleariella clade they range even 
higher, up to 27%. 
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